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Abstract: We present a new procedure, using on-shell recursion, to determine coefficients 
of integral functions appearing in the one-loop scattering amplitudes of gauge theories, in- 
cluding QCD. With this procedure, coefficients of integrals, including bubbles and triangles, 
can be determined without resorting to integration. We give criteria for avoiding spurious 
singularities and boundary terms that would invalidate the recursion. As an example where 
the criteria are satisfied, we obtain all cut-constructible contributions to the one-loop n-gluon 

scattering amplitude, .. [- + + .••), with split-helicity from an A/" = 1 chiral 

multiplet and from a complex scalar. Using the supersymmetric decomposition, these are 
ingredients in the construction of QCD amplitudes with the same helicities. This method 
requires prior knowledge of amplitudes with sufficiently large numbers of legs as input. In 
many cases, these are already known in compact forms from the unitarity method. 
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1. Introduction 

With the approach of the Large Hadron Collider, perturbative QCD computations will play 
an ever more prominent role in helping to unravel physics beyond the Standard Model. Com- 
putations at tree-level provide a ffist step for providing theoretical guidance. For reliable 
quantitative computations, at a minimum, next-to-leading order calculations are required. 
Over the past decade important strides have been taken in producing new next-to-leading 
order calculations to various standard model processes in pure QCD, heavy quark physics, 
and mixed QCD/electroweak processes. (See, for example, refs. [1] for recent summaries.) 
However, much more needs to be done to take full advantage of the forthcoming data. 
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In this paper we discuss one aspect of this problem: the computation of coefficients of 
integral functions appearing in one-loop scattering amplitudes with large numbers of external 
legs. In a very recent paper [2], an improved unitarity- factorisation bootstrap approach [3] 
was presented for obtaining complete amplitudes, including the rational function terms. As 
examples, one-loop six- and seven-point QCD amplitudes with two colour-adjacent negative 
helicities and the rest positive were computed. These rational function pieces are trivial in 
supersymmetric theories, but presented a bottleneck to further progress in QCD loop calcu- 
lations. With this obstruction mitigated, it is worthwhile to re-examine the cut containing 
pieces, as we do here, to see if further computational improvements may be found, especially 
in the light of recent progress. 

The recent advances have been stimulated by Witten's proposal of a 'weak-weak' duality 
between a string theory and J\f — 4 gauge theory [4] (generalising a previous description of the 
simplest gauge theory amplitudes by Nair [5]). (For a recent review of these developments, see 
ref. [6] .) Some of these efforts have been focused on analysing amplitudes from the viewpoint 
oiJ\f — A super- Yang- Mills theory as a topological string theory in twistor space [7]. There has 
at the same time been rapid progress for calculating amplitudes on the field theory side. This 
has led to a variety of calculations and insights into the structure of tree amplitudes [8-15] 
and supersymmetric one-loop amplitudes [16-32]. The novel organisation of string amplitudes 
has led to reorganisations of previously computed loop scattering amplitudes [33,34]. There 
have also been some promising steps for massive theories [35,36] and at multiloops [37]. 
The conclusion of these studies is that gauge theory helicity amplitudes are extremely simple, 
especially when compared to expectations. This simplicity has been the key to rapid progress. 

For the pieces of one-loop amplitudes containing cuts, the unitarity method [38, 39, 3, 
40,19,23,25], due to Dixon, Kosower and two of the authors, provides an effective means 
of computation. This method has been applied successfully in a variety of calculations in 
QCD [41,3,42,40] and in supersymmetric theories [38,39,37,19,25]. A recent improvement [23] 
to the unitarity method making use of generalised unitarity [3,40], allows coefficients of box 
integrals to be determined without explicit integration. Other related methods include the 
very beautiful application of MHV vertices [8] to loop calculations [16, 20] and the use [17, 
18, 21] of the holomorphic anomaly [43] to freeze [44, 17] the cut integrals. The unitarity 
method can also be used to determine complete non-supersymmetric amplitudes, including 
the rational function terms, by applying D-dimensional unitarity; a number of examples for 
four external particles, and up to six for special helicity configurations, have been worked out 
this way [41,45,42,40,32]. 

We make use of the decomposition of one-loop QCD n-gluon scattering amplitudes into 
contributions originating from J\f = A and A/" = 1 supersymmetric multiplets and those 
originating from a scalar loop [46,38,39]. Although supersymmetric multiplets contain a 
more complicated particle content, it is simpler to calculate the QCD amplitudes in this 
fashion; each of the pieces is naturally distinguished by their differing analytic structures. 
The simplest contribution is the most supersymmetric multiplet, i.e., the contribution of an 



- 2 - 



J\f — 4 super Yang-Mills multiplet. The Af = 4 amplitudes can be expressed as scalar box 
integral functions with rational coefficients [38] and these coefficients have been evaluated in a 
closed form for the case with maximally helicity violating (MHV) helicity configurations where 
two gluons are of negative hehcity and the rest of positive hehcity [38, 16] and for the case 
of next-to-MHV (NMHV) amplitudes [39,18,19,25]. (At tree-level the corresponding hehcity 
amplitudes are the Parke- Taylor amplitudes [47].) Using supersymmetric Ward identities for 
one-loop A/" = 4 NMHV box coefficients, expressions involving external gluino and quark 
legs have also been presented in ref. [31]. Results in the case of an A/" = 1 chiral multiplet 
have also been computed including, all-n one-loop MHV and all one-loop gluon six-point 
amplitudes [39,20-22,28,29]. In the scalar loop case the cut parts of the MHV amphtudes are 
known [38, 24] for an arbitrary number of legs. Very recently an explicit form of the complete 
amplitude, including rational functions has been presented at six-points for the case of nearest 
neighbouring negative helicities in the colour ordering [2]. 

In this paper we focus on the question of improving calculational methods for determining 
the coefficients of the integrals functions containing the cut constructible parts of amplitudes. 
As explained in ref. [23], by considering quadruple cuts in the unitarity method one can 
determine the coefficients of box integrals rather directly, without the needing to perform 
any integral reductions or integrations. However, in evaluations of A/" = 1 super-QCD or in 
non-supersymmetric QCD one encounters bubble or triangle tensor integrals, which have to 
be directly integrated; in general, experience shows that it is best to avoid direct integration, 
whenever possible. Even with this complication, as mentioned above, a wide variety of ampli- 
tudes in the J\f — 1 theory have already been computed. However, one may wonder if further 
improvements can be found. 

The approach we take here is based on recursion. Recursion relations have been very suc- 
cessfully used in QCD, starting with the Berends-Giele recursion relations [48]. The Berends- 
Giele recursion relations have also been applied in some special cases at loop level [49]. Re- 
cently at tree-level a new recursive approach using on-shell amplitudes, but with complex 
momenta has been written down by Britto, Cachazo and Feng [11]. This has led to compact 
expressions for tree-level gravity and gauge theory amplitudes [13-15]. This approach appears 
to work very generally, including massive amplitudes [36] . 

An elegant and simple proof of the on-shell recursion relations has been given by Britto, 
Cachazo, Feng and Witten [12]. Their proof is actually quite general, and applies to any 
rational function of the external spinors satisfying certain scaling and factorisation properties. 
Although the recursion relations emerged from loop calculations [19,50,11], the proof does 
not extend straightforwardly to loop level because of a number of issues. The most obvious 
difficulty is that loop amplitudes contain branch cuts which violate the assumption that the 
amplitude consists of a sum of simple poles, used to prove the recursion. Moreover, there are 
no theorems as yet on the factorisation properties of loop amplitudes with complex momenta. 
Indeed there are a set of 'unreal pole' [27,30] contributions that must be taken into account, 
but whose nature is not yet fully understood. 
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Initial steps in extending on-shell recursion to computations of loop amplitudes were taken 
in refs. [27,30], where compact expressions were found for all remaining loop amplitudes in 
QCD which are finite and composed purely of rational functions. This however, still left the 
question of how one might compute cut containing one-loop amplitudes in QCD via on-shell 
recursion. As mentioned above, a new approach has very recently been devised for doing 
so [2]. The method systematises an earlier unitarity-factorisation approach first applied to 
computing the one-loop amplitudes required for Z ^ 4 jets and pp ^ W -\- 2 jets at next-to- 
leading order in the QCD coupling [3]. With this method one assumes that the cut containing 
parts have already been computed using the unitarity or related methods. The new approach 
then provides a systematic means via on-shell recursion for obtaining the rational function 
pieces, accounting for overlaps with the cut containing terms. 

In this paper, we address the complementary question of whether one can also construct 
on-shell recursion relations for the cut containing pieces. We sidestep the issue of applying 
on-shell recursion in the presence of branch cuts by considering recursion for the rational 
coefficients of integral functions and not for full amphtudes. The key ingredients that will 
enter our construction are: 

• The universal factorisation properties that one-loop amplitudes must satisfy, allowing 
us to deduce the factorisation properties of integral coefficients [38,39,51]. 

• The structure of spurious singularities that can appear in dimensionally regularised 
one-loop amplitudes [52,51]. 

• The basis of integrals in which dimensionally regularised amplitudes may be expressed [53, 



• The proof of the tree-level on-shell recursion relations due to Britto, Cachazo, Feng and 
Witten [12]. 

• Prior calculations of integral coefficients that may be used as starting points in the 
recursive procedure. In this sense our procedure dovetails very nicely with unitarity 
methods, which provide the necessary inputs to the recursion. 

We will present simple sufficiency criteria to guarantee that a recursion on the coefficients 
is valid, with no boundary terms or spurious singularities affecting the recursion. 

As a non-trivial example of our recursion procedure we will consider one-loop n-gluon 
"spht-helicity" amplitudes of the form. 



where all negative helicity gluons are colour-adjacent. At tree level, amplitudes with this 
helicity configuration were computed in ref. [15]. In this paper we compute two components 
of the supersymmetric decomposition of one-loop QCD amplitudes for this configuration: The 



52]. 
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M = 1 chiral multiplet contribution and the non-rational or "cut-constructible" part of the 
scalar amplitude. For three negative helicities n-point amplitudes of the split-helicity form, 
the M = A [25] and M = 1 [28] amplitudes were derived earlier. We will present the cut parts 
of the scalar loop contributions, leaving the rational function parts as the final piece to be 
computed. 



2. Notation 

Throughout this paper we employ colour-ordered amplitudes. At tree level, the full amplitudes 
decomposes in the following way 

<^<'(l,2,...,n) = g^--'^ Tr(rMi)TM2)...TMn))ArKl),a(2),...,a(n)), (2.1) 

n 

where Sn/Zn is the set of all permutations, but with cyclic rotations removed. The T""^ are 
fundamental representation matrices for the Yang-Mills gauge group SU{Nc) , normalised so 
that Tr(T°T*) = 6°''' . (For more detail on the tree and one-loop colour ordering of gauge 
theory amplitudes see refs. [54,55].) 

For one-loop amplitudes of adjoint representation particles in the loop, one may perform 
a colour decomposition similar to the tree- level decomposition ( |2.1| ) [55]. In this case there 
are two traces over colour matrices and the result takes the form, 

Ln/2J+1 

{{h, a,}) = g"- E Gr„;e {a) , (2.2) 

C=l (TGSn/Sn;c 

where \_x\ is the largest integer less than or equal to x . The leading colour-structure factor, 

Gr„;i(l) = Nc Tr{T''' ■■■T"-) , (2.3) 

is just Nc times the tree colour factor, and the subleading colour structures (c > 1) are given 
by, 

Gr„;c(l) = Tr(T"^---T"^-^) Tr (T"= ■ ■ ■ T"") . (2.4) 

Sn is the set of all permutations of n objects and Sn-c is the subset leaving Gr^.^ invariant. Once 
again it is convenient to use U{Nc) matrices; the extra U{1) decouples [55]. The contributions 
with fundamental representation quarks can be obtained from the same partial amplitudes, 
except that sum runs only over the An-i and the overall factor of in Gr„.i is dropped. 

For one-loop amplitudes the subleading in colour amplitudes An-c , c > 1 , may be obtained 
from summations of permutations of the leading in colour amplitude [38], hence, we need only 
focus on the leading in colour amplitude An-i , which we will generally abbreviate to An , and 
use this relationship to generate the full amplitude if required. 
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We will represent amplitudes employing the spinor-helicity formalism [56]. Here ampli- 
tudes are expressed in terms of spinor inner-products 

(jO = {j-\l+) = u^{pj)u+ipi) , = (j + |r) = u+{pj)u4pi), (2.5) 

where u±{p) is a massless Weyl spinor with momentum p and positive or negative chirality. 
With the normahsations used here, [ij] — sign(p°p°) {j i)* so that, 

= '^Pi-Pj = Sij, (2.6) 

where 

= i^i+i = iPi + Pi+i? and ti^j = Kl- = (p. + . . . + p.f , (2.7) 

for generic cyclic ordered sets of external momenta j — pi + . . . + pj counting from leg i 
to end leg j . Note that [i j] defined in this way differs by an overall sign from the notation 
commonly used in twistor-space studies [4] . As in the twistor-space studies we define 

Xi = u+{pi) , Xi = u-{pi) . (2.8) 
We will also define spinor-strings such as 

[k\Ki,\l) ^ {k^\Ki,j\l^) ^ {l-\KiAk-) ^ mM ^ Y.[ka\{al) , (2.9) 

a=i 

as well as 

i n 

{k\K,^,Km,n\l) = {k-\f{ij^m,n\l^) = ^ (fc a) [a 6] (6 /) , (2.10) 

a=i b=m 

and 

j n 

[k\K,,,K^,n\l] = {k+\^^,,^m,n\n = • (2.11) 

a=i b=m 

In some cases, we also make use of commutators, such as 

{k\[q,K]\l) = {k\qK\l) - {k\Kq\l) . (2.12) 



3. Supersymmetric Decomposition of QCD Amplitudes 

The one-loop amplitudes for gluon scattering depend on the particle content of the theory 
since any particle of the theory with gauge charge may circulate in the loop. Let denote 
the contribution to gluon scattering due to an (adjoint representation) particle of spin J . The 
three choices we are primarily interested in are gluons (J = 1), adjoint fermions (J = 1/2) 
and adjoint scalars {J — 0). Instead of calculating contributions for these three particle types 
directly, it is considerably easier to instead calculate the contributions due to supersymmetric 
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matter multiplets together with that of a complex scalar (sometimes denoted as A/" = 
supersymmetry) . The three types of supersymmetric multiplet are the J\f = 4 multiplet and 
the J\f = 1 vector and matter multiplets. These contributions are related to the An^ by 



^Ar=l vector ^ ^[1] ^ ^[1/2]^ ^3^^) 

.Ar=lchiral ^ .[1/2] , .[0] 
— + ■ 

The contributions from these three multiplets are not independent but satisfy 

^A/'=l vector _ j^J\f=4 2^A/'=1 chiral (3 2) 

We can then invert these relationships to obtain the amplitudes for QCD via 

4 [1/2] ^ 4Ar=lchiral _ .[0] 



(3.3) 



Although the above relationship is for an adjoint fermion in the loop, the contribution from 
a massless fundamental representation quark loop can be obtained from the leading colour 
contribution of the adjoint case [46,39]. In this paper, we shall be focusing on two elements 
of this supersymmetric decomposition, namely the A^^'^ ^^^^'^^ and An^ terms. Together with 
the results for A^^^ amplitudes these form the components of QCD amphtudes. 

A one-loop amplitude may be expressed as a linear combination of dimensionally reg- 
ularised scalar integral functions Jj with rational (in the variables and A^) coefficients 
Ci, [53,52,38] 

A = Y,cJ,. (3.4) 

i 

In QCD, in general, after series expanding in e = (D — 4)/2, there are also finite rational 
function contributions arising from 0{e) terms striking ultraviolet poles in e. 

The above supersymmetric decomposition ( p.3|) of QCD amplitudes is useful because of 
the differing analytic properties of the components. Each of the components can be computed 
separately and then reassembled at the end to obtain QCD amplitudes [46,38,39]. For 
supersymmetric amplitudes the summation is over a smaller set of integral functions than in 
the generic non-supersymmetric set. In A/" = 4 theories, the functions that appear are only 
scalar box functions; whereas for A/" = 1 theories we are limited to scalar boxes, scalar triangles 
and scalar bubbles with no additional rational pieces. In ref. [38,39] it was demonstrated that 
supersymmetric amplitudes are "cut-constructible" using only four-dimensional cuts, i.e., the 
coefficients Cj can be entirely determined by knowledge of the four- dimensional cuts of the 
amplitude. The scalar loop contributions are, however, not fully constructible from their four- 
dimensional cuts and contains additional rational functions. Such terms arise, for example, 
from tensor bubble integrals. One may obtain these from cuts as well but the cuts must 
be evaluated in dimensional regularisation by working to higher-orders in the dimensional 
regularisation parameter e [41,45,42,40,32]. 
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4. Recursion Relations for One-Loop Integral Coefficients 



In this section we will develop tools for constructing the coefficients of the integral functions 
in one-loop amplitudes recursively. 

4.1 Review of Recursion Relations for Tree Scattering Amplitudes 

We first briefly review the proof [12] of tree-level on-shell recursion relations [11]. Let us begin 
by considering a generic tree amplitude A{pi,p2, . . . ,Pn) in complex momentum space and its 
behaviour under the following shift of the spinors, 

(4.1) 

Xb — > Xb — zXa , 

where z is an arbitrary complex variable. This will shift the momentum of the legs a and b , 

Paiz) = XaK + zXaXb: ^^2) 
Pb{z) = XbXb - zXaXb. 

By this shift, legs a and b remain on shell, Pa{z) = pl{z) = , and the combination pa ( 2;) +p;,( 2;) 
is independent of the parameter z . Under the shift, the amplitude becomes, 

^(Pl,P2,---,Pn) MPuP2,---,Pa(z),...,Pb(z),...,Pn) = A(z) , (4.3) 

where the shift respects the total momentum conservation of the amplitude. The shifted 
amplitude A{z) is an analytic continuation of the physical on-shell unshifted amplitude A{0) 
into the complex plane. 

For a tree amplitude in a gauge theory it has been proven in rcf. [12] that there will always 
be shifts for which the shifted amplitude A{z) will vanish as \z\ oo . Moreover, it can be 
shown that A{z) only has simple poles in z and the only poles which are present are the same 
ones which are present for real momenta. We can then evaluate the following contour integral 
at infinity, 

±.J^A{z) = Coo = ^(0) + 5]Res,=,„^, (4.4) 

a 

where Za are the simple poles of A{z) in z and ReSz=z„, signifies the residue at the location of 
the poles. Using the vanishing of A{z) as — > oo the boundary term Coo also vanishes and 

^(0) = -^Res,=,„^. (4.5) 

z 

a 

The poles and residues in the shifted amplitude A{z) are determined by the factorisation 
properties of the on-shell amplitude. An on-shell amplitude will factorise into a product of 
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two on-shell tree amplitudes as Kf- = (fcj + ■ ■ ■ + kjY — ^ , i.e. we will have the following 
factorisation 



A ~ — > A{ki, ■ ■ ■ , kj, K^j) X X A{kj+i, ■ ■ ■ , h^i, -Kij) . (4.6) 



At tree-level one can show that there are no other factorisations in the complex plane. Using 
these factorisation limits one can then write the following recursion relation for tree amplitudes 

^(0) = ^n-m„+l(^a)-^^mt+l(^Q) ' i^'^) 

a,h ° 

where A'f^_^^^^{za) and ^mt+i('^o) shifted amplitudes evaluated at the residue value Za , 
and h denotes the helicity of the intermediate state corresponding to the propagator term 



i/K'^ . In eq. ( |4.7D the sum selects tree amplitudes where the reference leg pa is contained in 



one tree amplitude and the reference leg pb is contained in the other. 
4.2 Recursion Relations for Integral Coefficients 

For the special set of one-loop amplitudes which consist of only rational functions (the "finite" 
loop amplitudes), a set of recursion relations were developed in ref. [27,30] for determining 
all such amplitudes. One-loop recursion relations amplitudes are more subtle than at tree- 
level because of the appearance of double and 'unreal' poles whose nature are not yet fully 
understood. More generally the presence of branch cuts invalidates a basic assumption of the 
tree-level proof outlined above. 

Very recently, a unitarity-factorisation bootstrap method was derived for dealing with the 
more general case containing cuts [2], systematising an earlier unitarity-factorisation bootstrap 
approach [3]. With this method the cut containing terms are evaluated using the unitarity 
method and the rational function pieces by finding an appropriate on-shell recursion relation, 
accounting for overlaps between the two types of terms. In practice many results are already 
known for the cut containing pieces, so in these cases only the rational function parts need 
to be determined in order to have complete QCD amplitudes. As an example, an explicit 
construction of the rational parts of the six-gluon QCD amplitude A^.^^(l~, 2^, 3+, 4+, S"*", 6~^) 
was given in ref. [2]. 

Here we wish to address the complementary question of whether we can also apply re- 
cursion relations to the cut containing pieces. The approach we take will be to construct 
recursion relations, not for amplitudes but for the coefficients of the integrals since these are 
cut free rational functions. In this way we avoid the issue of constructing recursion rela- 
tions in the presence of branch cuts. However, to be successful, a firm grasp of the analytic 
properties of the integral coefficients is required. These properties differ from those of full 
amplitudes: in general these coefficients will not contain the complete set of poles appearing 
in the amplitudes and they will contain spurious poles. The spurious poles can interfere with 
the construction of recursion relations since there are no theorems to guide factorisations on 
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spurious denominators. These spurious poles in general cancel only in the sum over all terms 
in the amplitude, and not in individual integral coefficients. But with an understanding of 
which spurious poles appear in which integral coefficients, we can maneuver around them. A 
detailed discussion of the factorisation properties of one-loop amplitudes, as well the spurious 
singularities that appear, may be found in refs. [38,51]. Here we will only summarise a few 
salient points. 

First consider the factorisation where a three vertex is isolated as one of the factors. For 
massless particles and real momenta such vertices vanish on-shell, though not for complex 
momenta [11]. The real momentum factorisation is more conventionally presented in terms 
of limits where the momenta of two legs Pa and ph become collinear [38,39,51,57,58], 



(4.8) 



+ V Split°7'°°P(a'*^ b''')A'^l\{. . . , (P)^ 



where P = pa + Pb ■ The case of multi-particle factorisation is also well understood; one-loop 
amplitudes factorise according to a universal formula [51], 



^one-loop i,i + m-l^ ^ ^ 



h=± 



/lone-loop/ jy^h \ ^ Atree f f jy- \-h \ 

I- • • ; ■:■■■) t^2 ^n-m+l\- • • ' I J^i,i+rn-lJ i ■ ■ ■) 

i,i+m—l 

_l_ ^tree / jv-Zi \ ^ <one-loop/ ( ^ IT \ 

"T ^m+l\- ■ ■ ? ^i,i+m-l ■>■■■) t^2 I- • • ) I ^i,i+m-lj ) ■ ■ ■) 

i,i+iri—l 

h ^ h '2 

+ ^m+l(- • • 5 ^i,i+m-l ; • • ■) ^n-m+l(- • • ? (~-^i,i+m-l) ; • • •) -^n(-^i,i+rri-li Pi' • • • y Pn) 

i,i-\-m—l 

{4.9) 

where -K^fj^i.i is the momentum invariant on which the amplitude factorises. The 'factori- 
sation function' denoted by jF„(Kj?j_,_„_^; pi, . . . actually represents a non-factorisation, 
since it contains kinematic invariants which cross the pole; they however have a universal 
structure linked to infrared divergences and are given explicitly in ref. [51]. 

The factorisation of amplitudes follows from the combined behaviour of integral functions 
and the integral coefficients in the factorisation limit. If we turn this around, given the general 
factorisation ( ^.9] ) of an amplitude and given its basis of integral functions, we may then 
determine the factorisation properties of the integral coefficients. 

Consider ffist the collinear limits. If we expand the amplitude in terms of integrals times 
coefficients using eq. (0), we have that 

E ^V./.n ^ Yl SplitlT(«'^% &'^^) E <n-lkn-l + ^ SplitT^°°^(«'% &'^)^n-T • 
i h I h 

(4.10) 

In order to disentangle the behaviour of the coefficients we need to know how the integrals 
flow into each other under collinear factorisation [38,51]. It is useful to choose a good basis of 
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integral functions with appropriate factors moved between coefficients and integrals in order 
to have a simple factorisation behaviour. We will make one such choice below when we discuss 
the examples. 



Now focus on the behaviour of a single term Ci^nh^n in the sum (|4.1CI| ) as two momenta 
become collinear. The very simplest case is when the two momenta pa and pb becoming 
coUinear belong to legs in a single cluster of external legs of the integral. In the case that 
the cluster has at least three legs so it does not become massless, the integrals have a smooth 
degeneration in the collinear limits as illustrated, 




a\\h 




where P = pa + Pb becomes massless in the factorisation limit. 
In this case the factorisation of the coefficient is simply 



Yl Split 



tree 
h 



i,n—l ; 



(4.11) 



where the coefficient on the left belongs to the left integral in the figure and the coefficient 
on the right belongs to the right integral. The coefficient has this simple behaviour because 
in this case there are no contributions to the Split°')^~'°°'' term in eq. ( 4.1(]| ). Contributions 
from this term all come either from discontinuities as massive clusters degenerate to massless 
ones or from integrals with two massless legs carrying the collinear momenta pa and pb ■ (See 
ref. [51] for further details.) 



By applying the same logic to multi-particle factorisations ( |4.9D we conclude that the 
coefficients also behave as if they were tree amplitudes as long as the factorisations are entirely 
within a cluster of legs (and are not on the momentum invariant of the entire cluster). That 
is the coefficient behaves as, 



h 

;,m+l ■ 



(4.12) 



Assuming that the spurious denominators do not pick up a ^ dependence — below we 
describe simple criteria for ensuring this — we obtain a recursion relation for the coefficients 
which strikingly is no more complicated than for tree amplitudes. 



Cn(0) = 



a,h 



^2 '^m.a+l\^o 



(4.13) 
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where and c^_m„+i(-2^a) shifted tree amphtudes and coefficients evaluated at 

the residue value , h denotes the helicity of the intermediate state corresponding to the 
propagator term i / . In this expression one should only sum over a limited set of poles; if 
the shifts are chosen from within a cluster, the only poles that should be included are from 
within the kinematic invariants formed from the momenta making up the cluster. Pictorially, 
this coefficient recursion relation is, 




In order to have a valid bootstrap we must have that the shifted coefficient vanishes as 
\z\ oo; otherwise there would be a dropped boundary term. We can, however, impose 
criteria to prevent this from happening. Consider an integral and consider the unitarity cut 
which isolates the cluster on which the recursion will be performed, i. e. the one with the two 
shifted legs. 




The dashed line in this figure indicates the cut. The recursion is to be performed with the 
two shifted legs from the right-most cluster. Then simple criteria for a valid recursion are: 

1. The shifted tree amplitude, on the side of the cluster undergoing recursion, vanishes as 
\z\ — >• oo . 

2. All loop momentum dependent kinematic poles are unmodified by the shift ( [4.1|) . 

Since the location of none of the poles in z depends on loop momentum and any shifted 
tree amplitude can be written in the form of poles in z times residues (which are of course 
independent of z), the z dependence appears only as an overall prefactor in front of the 
loop integral. We are therefore guaranteed that the shifted coefficient will have vanishing 
behaviour as \z\ oo . With these two criteria we have an added bonus: The spurious 
singularities appearing in the shifted coefficients cannot depend on z since the integration 
does not tangle with z . This guarantees that our recursion relations are not contaminated 
by spurious singularities. In a direct integration of the cuts, spurious singularities would 
arise from evaluating the integrals, but since z-dependence does not enter into the integrals 
they cannot enter into spurious singularities. All calculations that we perform in this paper 
satisfy these two conditions proving that there are no boundary terms or spurious singularities 
invalidating the recursion. 
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The above criteria are simple sufficiency conditions for constructing a recursion for loop 
integrals, but are not necessarily required: It is possible to find valid shifts outside the above 
class, for example, involving shifted legs from different corners of triangle or box integrals, 
but we will not discuss this here. In general, one would also need to account for factorisations 
involving the second term in eq. ([4.10|) as well as non-trivial flows under factorisation of 



integral functions into each other and into loop splitting and factorisation functions [51]. 
In addition some care is required to avoid shifting the spurious denominators. We leave a 
discussion of these issues to future studies. 

4.3 Example 

To make our procedure more concrete, it is useful to illustrate it with an explicit example, 
before turning to the main calculation of this paper. In this section we will use recursion 
relations to compute the coefficients of the six-point split-helicity one-loop amplitudes in the 
= 1 chiral multiplet as well as for the scalar loop contribution A^^^}-^ . 

In general, one-loop amplitudes contain a set of integral functions. For the amplitudes of 
interest here, a convenient set is [46], 



Ko(r) 



e l-2e 



(^-log(-r) + 2 + ^) + 0{e) 



log(r) T ^ ^ - ^^Sl'^) - (r - l/r)/2 



(4.14) 



1 — r (1 — r 



|3 



We find it convenient to use the Lq , L2 and Kq functions to build up the expressions for the 
complete amplitudes. This choice of integral functions is a good one from the viewpoint of 
having simple factorisation properties for the integral coefficients. The Lq and L2 functions 
can be thought of either as the difference of two bubble integrals or as a two-mass triangle with 
an integrand which is linear in the Feynman parameter for the leg linking the two massive legs. 
The Lj functions have the advantage of being finite in e and free from spurious singularities 
as r — i> 1 . 

The dimensionally regularised amplitudes carry an overall factor of cr , 

1 r(i + e)r^(i-e) 

'' = J4^ r(l-2e) ' 

which we suppress in this paper. 

In our example, we will use a five-gluon one-loop amplitude with an A/" = 1 chiral multiplet 
in the loop as the starting point in the recursion. This five-gluon amplitude is [46], 

^Af=l chiral^^-^ 4+, 5+) = ^ A'^^%1-, 2-, 3-, 4+, 5+) ( Ko(t5,l) + Ko(t3,4)) 

+ i c^=^ "(5+, 1-; 2-; 3^, 4+) M^Mt^!:!!! , 



2 



''5,2 

(4.16) 
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where, 

^tree.^- 2- 3- 4+ 5 



[4 5]= 



[51] [12] [2 3] [3 4]' 

^ ' ' ' ' ^ [51] [12] [2 3] [3 4] 

In general we will use semicolons to delineate the groupings of the legs into clusters of the asso- 
ciated triangle integrals. In this case the triangle diagram representation of Lo(— t5,i/(— ^5,2)) 
is 




For the cut-containing terms of the scalar loop contribution to the five-gluon amplitude 



A[°l(l-,2-,3-,4+,5+) = -A-^=i"(l-,2-,3-,4+,5- 

3 



+ \ cM(5+, 1-; 2-; 3", 4+) ^ ~ ,'!'^//3~'^'^^^ + rational . 
3 (^5,2)-^ 



(4.18) 



where the coefficient c'"' is given by, 

^ ' ' ' ' ^ [51] [12] [2 3] [3 4] ^ ' 

In this case the rational terms in the amplitude are known [46], but we will not need these 
here. (A very recent discussion of rational parts is given in ref. [2].) Note that the two 
coefficients c-'^^^ ^^"^^ and c'''' have different dimensions, because they multiply factors of 
different dimensions. (The functions Lq and L2 are dimensionless.) 

The example here illustrates how recursion relations for loop coefficients work and will 
lay the ground for the generic computation of the n-point split-helicity coefficients c^^^ ^^"^^ 
and Cn' . As we will see, the coefficients c:lf^^ ^^"'^^ and Cn^ are computed as sums over different 
shifts of the five-point coefficients c^^^ '^^"'^^ and cf' . It will turn out that these coefficients 
are quite similar and we can therefore unify the expression for the coefficients using a variable 
^ and ^ which both takes the identity value I for the A/" = 1 chiral coefficients, where I is 
defined to be the identity matrix in the space of gamma matrices and Jf = ^2-^5,1 , = -^5,1 ^2 
in the scalar loop contribution. The five-point 'googly' coefficients C5 for both theories can 
then be written as, 

c. - cf5+ 1-- 2 - 3 4+) - J4|^|5][4|P|5][4|[fc2,i^d|5] 20) 
C5-c(5 ,1 ,2 ,3 ,4 ) - ^ [5 1] [1 2] [2 3] [3 4] " ^ ' 
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(In twistor terminology, 'googly' means the conjugate amplitude where positive and negative 
helicities are interchanged. In the MHV case, a googly amplitude has two positive and the 
rest negative helicities.) 

Now consider the six-point coefficient cq = cq{6^, 1^; 2^; 3^, 4+, 5+) of the integral func- 
tions Lo(— t6,i/(~^6,2)) and L2(— t6,i/(— ^6,2)) • Graphically the associated integral function 
is, 



6 




In this example we will use the shift. 



— > A'? — zXa , 

(4.21) 

A4 — > A4 + zAs . 

We will first check that the two criteria for finding valid shifts are satisfied. Since the shifted 
legs are from the {3, 4, 5} cluster, we must check the cut 

2- 




6+ 5- 

For the case of a scalar circulating in the loop the tree amplitude isolated on the right side of 
the cut is. 



where the subscript s denotes a scalar circulating in the loop. Under the shift ( [4.21| ), the tree 
amplitude becomes, 

Atree(o- 4+ c:+ p+ f > - . y\ _ ,/ (3 4) (3 ^j) ( A 

^ ' ' '^'^'^ (3 4)((4 5) + 2(3 5))(5 4)(4£i)(£i3)' ^^''^^ 

and we immediately see that both criteria are satisfied: the z-dependence factors out of the 
integrand and the shifted coefficient times the integral vanishes as \z\ — > oo , because an 
overall prefactor vanishes. Similarly, it is easy to check that the criteria are satisfied also for 
a fermion circulating in the loop. (The M = 1 chiral case is just the sum of the fermion and 
scalar loop contributions so it too satisfies the criteria.) 

From eq. ( [4.23| ) it is clear that the only channel which can contribute to the recursion is 
the (4 5) channel, which is determined by the coUinear factorisation, 

C6(6+,l-;2-;3-,4+,5+) Splitl^^^(4+, 5+) 05(6+, 1"; 2"; 3", (4 + 5)+) . (4.24) 
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Thus the coefficient Cq{z) = c{pQ,pi;p2;p3{z),p4{z),p5) has only one pole in z. Notice that 
the coefficient 05(6+, 1^; 2^; 3^, (4 + 5)+) is given by the five-point amplitude ( [4.20| ). 

Following the tree-level construction described in ref. [11], the coefficient Cg can hence be 
computed as the residue of —cq{z)/z at the coUinear pole zi = — (4 5)/(3 5) , 



ce = cq{0) = —Res 



Z = Z\ 



Applying the recursion ( ^.131 ), 

= -(4 5)/(3 5), 
[4i^4,5] = [4|ir4,5|3)/cu, 

[2 3] = [2 3]-z[2 4] = [2|ir3,4|5)/(3 5) 

the triangle coefficient evaluates to, 



UJUJ 
[5i^4,5] 
[3it4,5] 



(3|i^4,5|4] , 

[5|i^4,5|3)/cu 



(4.25) 



(4.26) 



C6 



c(6+, I'- T- 3-, K+) — A(4+, 5+, (-ir45) 

■S45 



[ii^45|P|6] [i^45|P|6] 


i^45|[A;2,i^6,2]|6 


% (-0[4 5]3 


[61][12][2 3][3i^45] 


S45 [(-i^45)4][5(-i^45)] 



(3|i^3.5P|6] (3|ir3,5P|6] (3|jr3.5[fc2,i^6.. 

[2|ir3,5|5) [6 1] [12] (3 4) (4 5) ^3,5 



(4.27) 



The above expression for the six-point coefficient for M = \ agrees with the result first given 
in ref. [21], after setting Jj* and If equal to unit matrices. In the scalar loop case. If = ^2^q,i 
and f = ^6,1 h in eq. (14271 ). 



5. All-n Split-Helicity Integral Coefficients 

In this section, we present the coefficients of all integral functions in the split-helicity ampli- 
tudes, 

A(l-,2-,- ■■,/-,/ + l+,---n+), (5.1) 

with either an = 1 super-multiplet or scalar circulating in the loop. Within the super- 
symmetric decomposition of QCD amplitudes ( p.3| ), the TV = 4 component always consists 
entirely of box integrals [38] . For the split helicity configuration it is not difficult to show that 
for an A/" = 1 chiral multiplet or a scalar in the loop there are no boxes and the amplitudes 
consist entirely of triangle and bubble integrals [21]. 

The starting expressions in our recursion are the one-loop five-point split-MHV amplitudes 
which were first obtained in ref. [46]. As before, we use the Lj and Kq as the basis functions. 
From an examination of the unitarity cuts it is not difficult to show that the only logarithms 
that appear in the amplitude, to finite order in e, have arguments ti,j where 1 < i < I and 
I + 1 < j <n: otherwise a vanishing tree amplitude appears in the cuts. 
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Before continuing, it is convenient to employ a unified notation for both the M = 1 and 
scalar loop cases, so that we can present the computations in the two cases simultaneously. 
We define, 

P = 1, P = I, for the chiral A/" = 1 loop contribution , 

^ = fml/^L , = I/^L , for the scalar loop contribution , 

where I is the identity matrix, km is the momentum of the massless corner, and Kl denotes 
the total momentum of the left massive corner of the triangle. 




We write the amplitude for an A/" = 1 chiral multiplet in a Lq and Kq basis so that, 
A^=' "(r, 2-,- ■■/-,/+ 1+,- ■■ = aiKo(s„i) + a^K^isu+i) 



1 X ^ m,r Lq [tm,r/^m+l,r] 1 V ^ -r?i,r 



Lo[^m,f/^m+l,f] (5-3) 



9 / J n,L J. 9 / J n,i j. _ 

^ m,r Wl,r ^ Im+l.f 



where ^ = ^ = I in the Af = 1 triangle coefficients c^f and c^f . For the scalar loop 
contribution, the Kq functions and Lq functions are just equal to one-third of the J\f = 1 
contributions, as can be easily demonstrated from our recursion, so that we write the full 
amplitudes as 



■ n 



2-, ■■■,/-,/ + 1+, ■■ =1 A^=\l-, 2-, ■ ■ ■ , r , / + 1+ 

1 m,r '^2[tm,r/tm+l,r] 1 _m,f L2 [tm,f /^7n+l,f ] . , i^-^) 

+ oZ^ '^n,i 13 ^ qZ^ +3 ^ rational, 

where now ^ = ftm^L , ^ = for the triangle coefficients c^f and c™;'^ . In these 

expressions, c^f is the coefficient of the triangle function. 




while the c^f is the coefficient of the triangle function. 
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(to-1)+/\ (?fl + l)+ 

il + 1)^A /o/l\ A 

I- 

(f+l)~ f- 

Note that the coefficients (x^'l are simply related to the coefficients c™f by parity and rela- 
belling. Thus, we need only obtain the c^f coefficients to have a solution. 

The coefficients of the Kq functions satisfy the same recursion relations as the tree ampli- 
tudes with the starting point being a five-point tree amplitude. This gives us the coefficients 
of the Kq functions to be proportional to tree amplitudes, 

ai = a2 = \a''''\\-. ...,/-,/ + 1+, n+) , (5.5) 

The formulae for the tree amplitudes may be found in ref. [15]. 

As a warm up before dealing with the all-n point cases of triangle coefficients, we ffist 
consider specific and detailed examples of how to recursively derive the coefficients of the 
split-helicity triangle integral functions. 



5.1 Seven Points 

The ffist example we consider is the googly NMHV amplitude A(l~, 2~, 3~, 4~, 5^, G"*", T"*") 
and the coefficient c-j = 07(7^, 1~; 2~; 3~, 4~, 5+, 6+) which is associated with the integral 
functions Lo(-t7,i/(-t7,2)) and L2(-t7,i/(-t7,2)) , 

2- 

.3-_ 
6+ 

In the recursion we shift = p^^z) , p^ = p^i^z) according to 

A4 — > X4 — zXk, , 

(5.6) 

A5 — ^ A5 + 2A4 . 

The momenta p1{z) and pf{z) remain on-shell, i.e. p\{z) = and pl{z) = and momentum 
conservation is unchanged. The shift we consider here satisfies the two criteria for valid shifts 
given in section ^. We will demonstrate this for all-ra in section |5.3| . 

Under this shift, the contributing factorisations are given by the collinear limits for p-^ \\ 

C7(7+, 1-; 2-; 3", 4", 5+, 6+) Split+(3-, 4") 0^(7+, 1"; 2"; (3 + 4)", 5+, 6+) , (5.7) 
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and for \\ Pe , 

C7(7+,l-;2-;3-,4-,5+,6+) Split_(5+, 6+) C6(7+, 1"; 2"; 3", 4", (5 + 6)+) . (5.8) 

Thus the coefficient C7{z) has only two simple poles in z , corresponding to the two collinear 
factorisations. The six-point coefficients that go into the collinear factorisations p^ || p4 and 
P5 II Pe can both be derived by recursion from the five-point coefficient as we have seen in 



section 4.2 



By recursion we then calculate the coefficient C7 as the residue of —ci{z)/z at the poles 
zi = [3 4]/[3 5] and Z2 = -(5 6)/(4 6) , 

C7 = C7(0) = - V Res,=,„^, (5.9) 

z 

a=l,2 

i.e. the contribution from Zi is given by 

crU = C6(7+, 1-; 2"; K,^^ , 5+, 6+) — ^(3", 4", {-Ks,,)-) . (5.10) 

-534 

When inserting the form eq. (f4.27|) for Cq{7^, 1^; 2^; K^4^, 5^, 6"*^) we can write double 
shift of the five-point googly coefficient C5 , 

CjU, = C5(7+,l-;2-;4":4,^s"6)— ^(^+,6+,(-%)-)— A(3-,4-,(-i^3,4)+). (5.11) 

The expression for cyl^j is thus given by a sequence of shifts of the five-point googly coefficient 
C5 . The ffist shift adds a positive helicity leg, while the second shift adds a negative helicity 
leg. As we increase the number of legs the repeated shifts will be very useful as a way of 
organising the results. In the following we will call shifts that add a positive helicity leg a 
'plus shift' and a shift that adds a negative helicity leg a 'minus shift'. The double hat in the 
above equation expresses that we are shifting a leg twice. 

Explicitly we find for cjIz-^ , 



I ^ . [517^3,5^^3,6^17] [517^3,5^^3,6^17] [5 1 7^3,5^^3,6 [^2, ^^l] 1 7] . . 

""''^^ ' [3|K3,5|6)(6|ir3,6|2][71][12][34][45]t3,5t3,6 ' ^' ' 

The contribution Ci\z2 is given by a minus shift followed by a plus shift of C5 



crU = C5(7+,l-;2-;i^3-.,4t6)— ^(3~,4',(-i^3,4)^)— ^(5+,6+,(-i:5,6) 

' -534 556 

^ . (417^4,6^17] (4|ir4,6P|7] {A\K^,e[k2,KL]\7] 
' [3|i^4,6|6) [71] [12] [2 3] (4 5) (5 6) t4,6 



(5.13) 



We find it useful to keep track of the various contributions to the integral coefficients via 
representations of these as different paths in a helicity diagram. Such a diagram provides a 
pictorial expression for the factorisation structure of a massive corner in a coefficient. 
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In our conventions the point (1, 1) will correspond to a corner with one positive helicity 
leg and one negative, i.e. {p,n) = (1, 1) . The point (1,1) is thus associated with the googly 
five-point coefficient C5 . The path going one step up and right (/") corresponds to a plus 
shift, while a path going one up and left (\) should be associated with a minus shift. Hence 
the coefficient C5 is related by a plus shift to the part given by the googly six-point coefficient 
represented by (1,2). By a subsequent minus shift one can then get the part C7I21 of the 
seven-point coefficient corresponding to the point (2, 2) in the diagram. Another possibility 
is to go from the coefficient C5 to the NMHV six-point coefficient represented by (2, 1) by a 
minus shift and subsequently by a plus shift to get the part cyj^^ . 

The helicity diagrams corresponding to these two different shift paths are depicted below. 

(2,2) (2,2) 

(1,1) (1,1) 



C7U2 C7I 



A diagrammatic approach was also used at tree level to represent different contributions to 
split-helicity amplitudes in [15]. 

We conclude this subsection by noting that the complete expression for the googly NMHV 
coefficient C7 is given by the sum over shifts of the coefficient C5 . Expressed in terms of the 
helicity diagram, the coefficient is given as a sum over paths. We have checked that the 
expression for the seven-point-googly NMHV coefficient in the J\f = 1 case 

C7 = C7|^, + C7|^2 , (5.14) 

agrees with the results given previously in ref. [28]. 
5.2 Eight Points 

It is also worthwhile to work through an eight-point googly NMHV triangle coefficient since 
this case has some non-trivial aspects in its recursive build up which will be important for us 
when considering the generic ra-point triangle coefficients. 




In this case we use the shift 



A4 — ^ A4 — zX^ , 

A5 — A5 + 2;A4 . 
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Again, this shift satisfies our criteria for vahd shifts, as we will show below in section E]3 
The eight-point coefficient Cg = c(8+, 1^; 2^; 3^, 4^, 5+, 6+, 7+) is given by the sum 

Cg = Cglzi + C8I22 5 

csl,, = c(8+,l-;2-;i^3-4,5+,6+,7+)^A(3-,4-,(-i:3,4)+), (5.15) 

csl., = c(8+,l-;2-;3-,4-,K5+6,7+)— A(5+,6+,(-i^5,6)"). 

When inserting the expression for the seven-point coefficients into the recursion for cg , 
one finds that it is related to the coefficient C5 by three different sequences of shifts: 1) two 
plus shifts and a subsequent minus shift, 2) a plus shift, a minus shift and a plus shift and 
3) by a minus shift and subsequently two plus shifts. Schematically this means that the 
eight-point coefficient is given by 

Cg = cl + cl + 4 = C5(+ +) X A{+ + -) X A{+ + -) X A{- - +) 

+ C5(+ +) X A(+ + -) X A{- - +) X A{+ + -) (5.17) 

+ C5(+ +) X A{ +) X A{+ + -) X A{+ + -) . 






Co 



Again we can picture the various contributions to the eight-point coefficient in helicity dia- 
grams. 

5.3 All-n Triangle Coefficient 

We now discuss the recursive calculation of the n-point triangle coefficients c^'[ for the cases 
of a chiral Af = 1 multiplet or scalar in the loop. As before we calculate both of these cases 
simultaneously, using the unified notation of having P and P , defined in eq. (^]2|), to keep 
track of the two different types of loop contributions. 

A generic n-point triangle function in the split-helicity amplitude is diagrammatically. 




(r + 1) 
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At least one of the three-legs must be massless because otherwise the triple cut would vanish. 
(The case where the massless leg is of positive helicity can be obtained by conjugation.) 

Consider first a recursion on the right-most corner of the triangle. We shift legs / and 

This corresponds the shifts considered in our previous examples. The momenta pf{z) and 
pf+iiz) remain on-shell, i.e. pf{z) = and pf_^^{z) = and the total momentum remains zero. 



To check if this shift is valid according to the criteria given in section [4.2| , we need to 
consider the behaviour of the tree amplitudes isolated by the cut, 




(r + iy 

The tree amplitude on the left-hand side of the cut is, 

^'"'((-^//s)^, {r + 1)+, . . . , n+, 1-, . . . , m-, (-f^/J) , (5.19) 
and on the right-hand side is, 

A*-^(£;%, (m + 1)-, . . . , r , (/ + 1)+, . . . , r+, ijp , (5.20) 

The subscripts f/s signifies either a scalar or a fermion in the loop. 

Note that the shift does not modify the loop momenta appearing in the cut; if we take 
£ as the independent variable, i ' will not depend on z , as it is a function of the sum of the 
momenta in the tree, 

-i' = pm+i + ...+ Pi{z) + pi+i{z) + ...+ i, (5.21) 



which is independent of z because of eq. (|4.1| ). Thus the z-dependence cannot enter into the 
tree amphtude ( ^.19|) on the left-hand side of the cut. 

Now consider the tree-amplitude (|5.20| ) on the right-hand-side of the cut, which is the one 
containing the shifted legs. The shifted amplitude can be written as a sum of terms where 
each contains only a single pole in z corresponding to a factorisation of the tree amplitude. 
Only factorisations where loop-momentum flows between the two factor tree amplitudes can 
introduce a problematic pole, 

{i + ...+Pi{zW ^^-^^^ 
However there is no such factorisation of the tree, as the only potential factorisation that 
could contribute, 
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/T^ / + 1 

vanishes because the hehcity structure forces at least one of the two tree amphtudes appearing 
in the factorization to vanish, assuming there are at least a total of three external legs (not 
counting the ones carrying loop momentum). The vanishing of this pole contribution is 
related to the lack of box integrals in these amplitudes. This shows that the second criterion 
of section that all loop momentum dependent poles are unmodified by the shift is satisfied. 



The first criterion for the shift that the coefficients c(z) vanish when z goes to infinity 
is also satisfied. This follows from the vanishing of the tree-amplitude in the cut as z goes 
to infinity and because all z-dependence can be pulled out of the loop momentum integral. 
The limit \z\ —>■ oo then commutes trivially with the integration and we conclude that the 
coefficient times the integral vanishes in this limit. This proves that our recursion is valid and 
that there are no boundary terms. 

It is worth noting that since the shift comes out of the integration, all logarithms will be 
unaffected under the shifts, so that \n(tij) In(tjj). 

Since we have proven that our shift is valid for all integral coefficients in the split-helicity 
amplitudes, we may now consider the general recursion. With our choice of shift, the only 
relevant factorisations are given by the collinear limits for || pi , 

i-i\\i 



and for pi+i \\ pi+2 ■ 



m,r 



m,r 



l+l\\l+2 



Splits ((/-i)-,r)xCVi' (^-23) 



Split_ ((/ + !)+,(/ + 2)+) xC,,. (5.24) 



The coefficient c^f can now, as in the examples, be computed as the residues of —c^f{z) / z 
at its two collinear poles Zi = [I — 1, /]/[/ — 1, / + 1] and 2:2 = — (/ + 1, / + 2)/(/, / + 2) , 



c, 



= CfiO) = -J2 Res.=,^^ , (5.25) 

i=l,2 



which gives 




^n,l ~ 


m,r 1 
^n,l Ul 


•m,r 1 


m,r 


m,r 1 


m,r ( 
^n-l,l \-- 



I m.r 



nl \Z2 1 



I 



^1+1,1+2 



Sl-l,l 

A((TTT)'',(/ + 2)+,(-i^,+i,,+2)-). 



(5.26) 
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Thus, we have succeeded in expressing n-point coefficients in terms of (n— l)-point coefficients. 
The analogous steps can be apphed to the left massive corner. By repeated apphcation of the 
above steps c^f can be written in terms of shifts of the basic coefficient C5 . On each recursion 
the coefficient is a sum of two terms: one based on the {n — l)-point coefficient with one less 
negative helicity and the other based upon the n — 1 point coefficient with one less positive 
helicity. One can think of the recursion as leading via a "path" back to the coefficient of the 
five-point amplitude where the massive corner contains exactly one negative and one positive 
helicity leg. The coefficients c™f are then given by the sum, 



(5.27) 



The Tpj^ p^ term are generated from the same starting point; the five-point googly triangle 
coefficient c(5+, 1~; 2~; 3~, 4+) by sequences of plus/minus shifts on the left and right massive 
corners independently. The sequence of shifts is encoded in the arguments Pl and Pr and 
correspond to a solution of a set of parameters a^, pi, Pi', ai' in a helicity diagram. Two 
examples of paths for a massive corner are. 



ON 




PN+1 



UN 




Both paths contribute to the sum in eq. ( |5.27D of a coefficient with x negative helicity legs 
and y positive helicity legs on one corner. The coefficient is represented by the point [x, y) , 
with {x,y) = {I ~ m,r — I) . Since each massive leg must be reduced recursively, one sums 
over two sets of paths; the paths in a helicity diagram for the left corner and the paths in a 
helicity diagram for the right corner. Following this logic we have obtained a general solution 
for the split helicity with either an A/" = 1 multiplet or scalar in the loop. A path in a helicity 
diagram can be expressed as. 



Pl 
Pr 



PWhPi', 



N' + k' 



(5.28) 



where k, and A^, A^' are determined by the path in the helicity diagram. The parameters 
K, k' will either take the value or 1 depending on if the path in the helicity diagram starts 
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to the left or to the right. The parameters and A^' are given roughly by the number of 
corners of the path and can take values between 1 and the number of negative and positive 
legs whichever is least. In this range only those values for N and N' will contribute for which 
the inequalities eq. (|5.3CI| ) have solutions. 



The boundary parameters for specific k, k' and A^, A^' are 

= m + 1, On = I, Pn+k = 1 + 1, pi = r , 

(5.29) 

P^, = m - 1 , Pjy, = aN'+K' = n, (Xi = r + 1 . 

For a fixed set of boundary parameters the remaining parameters will have to satisfy the 
following set of inequalities, 



K, = 0,1, Oii < Oj+i, i = N — 2 + K , 

Pi > Pi+i, i = 1, A^ - 1 



5 



k' = 0,1, /3, > J = k',...,N'-2 + k', 

CTj < aj+i, j = 1, A^' - 1 . 

For a solution to the inequalities the Tp^^p^ takes the following compact form 

= ^(-1) 



(5.30) 



i^n+N' {mRPQL\l) {1\QrPQl\1) {l\QR[K^,K]QL\l) (r,r + 1) 



[1 2][2 3]...[/ - 1, /](/, / + l)...{n - 1, n){n, 1) 
nil [«i - 1, tti] (a, A + 1) 



X 



{pN\KN\aN - 1] nj=i {Pi\Ki\ai - 1] {ai\Ki\piJ^i + 1) 



1 

X 



(5.31) 



where 





= K 




= K 


Kr 


— K 


Qr 




P 


= kmKR,l, 









= -^/3j,crj + i , 






Ql 


— K' K' K' 


P 


= Kpkm,!. 



(5.32) 



In order to calculate a triangle coefficient one uses the above formula and sum overs the 
contributions to Tp^^p^ for each solution of the parameters ctj, Pj' and ai, pi'. This is equivalent 
to summing over all possible paths in a helicity diagram. 
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6. Results for NMHV Split-Helicity Amplitudes 



In this section we specialise the results of the previous section to the case where there are 
exactly three nearest neighbouring negative helicities in the colour ordering. For this "next- 
to-MHV" (NMHV) helicity configuration almost all of the super symmetric decomposition 
of the QCD amplitude is known: the — A contributions are known from ref. [25] as are 
the contributions from the M — 1 chiral multiplet [28]. The cut containing scalar loop 
contributions are contained in the expressions of the previous section and in the general 
formula for Tp^ p^. This leaves only the rational parts as undetermined. With an evaluation 
of these rational functions parts using, for example, the method of ref. [2] we would have the 
complete QCD amplitudes for this helicity configuration. 

The expressions for the M = 1 and scalar loop contributions are rather similar, whereas 
that of the A/" = 4 multiplet is very different: the A/" = 4 amplitude consists entirely of box 
functions. In contrast, the J\f = 1 contribution consists entirely of the Lq and Kq functions 
while the scalar loop contribution also contains L2 functions. 



6.1 Explicit Results 

Using the generic formula for Tp^ we have obtained the following compact results for the 
NMHV helicity configurations: For the amplitudes with an A/" = 1 chiral multiplet running 
in the loop the result is. 



^M=i chi.ai(i-^ 2-, 3- 4+, 5+, • • . , n+) = ^ (Ko(.„0 + Ko(.34)) - ^ £ l,r ^^^^flM 

r=A 2''' 

_ i_ ^ Lo[^2,r/^2,T-+l] _ _^ ? Lo[^3,r/^3,r+l] 

^ ^ '^2,r+l ^ ^ I3,r+1 

(6.1) 

where. 



dn,r 



{3\Ks,rK2,r\ 


|l)'(3|ir3,r[A:2,i^2,r]i^2,r| 


ll) 


[2|i^2,.| 


k)[2|^2,r| 


r + 1) (3 4) . . . (r - 1 r) (r + 1 r + 


2) . . . (nl) t2,rt3,r 



_ ^ (3|K3,,+3^2.,+3|l)' (3|i^3,,+3i^2a+3 [K+I, i^2,r] |1) (J + 3 J + 4) (5.2) 
gn,r — 



- [2|ir2,,+3|j + 3)[2|i^2,i+3U + 4) (34)(45)...(nl) tsj+st2. 



,i+3 



hn,r — ( l)"^n,n-r-+2|(i23..n)-»(321n..4) ■ 



This expression for the M —1 chiral amplitudes agrees with the one first given in ref. [28] . 
The expression for the cut-constructible parts of the scalar amplitude have a very similar 



- 26 - 



form apart from the appearance of L2 functions, 

2-, 3-, 4+, 5+, ■ ■ ■ , 72+) = 1 A^=' 2-, 3-, 4+, 5+, ■ ■ ■ , n+) 

_ i 7 L2[t3,r/i2,r] ^ ^ L2 [t2,r/i^2,r+l] ^ ? L2 [t3,r/i!^3,r+l] 

"J ^.r "J ^2,r+l ^=4 ''3,r+l 

+ rational , 

(6.3) 

where in this case, 

' _ (3|i^3,rfc2|l)(3|fc2i^2,r|l)(3|ir3,.[fc2,ir2,.]ir2,,|l) 



[2|if2,r|r)[2|ir2,r-|r + 1) (3 4) . . . (r - 1 r) (r + 1 r + 2) . . . (n 1) ' 

(3|ir3,,+3^2,,+3P|l) {3\Ks,j+sK2,j+3P\l) {3\Ks,j+3K2j+3 [K+i, K^^r] |1) (j + 3 J + 4) 



5'n, 



[2|i^2,i+3|j + 3)[2|ir2,,+3|j + 4) (34)(45)...(nl) t3,,+3t2,,+3 



hn,r -( 1) fi'ri,n-r+2|(i23..n)^(321n..4) ' 

(6.4) 

and P = fcr+i-ft'r+i,! and P = Kr+i^iK+i . 

Note that for the M = 1 case the infrared-singularities lie entirely within the Kq functions 
with the expected result, 

A-^=i(l-,2",3-,4+,5+,---,n+)|,-i = ^xA*=^^^(l-,2-,3^4+,5+,---,n+). (6.5) 

6.2 Obtaining the NMHV Coefficients from General Formula 

In this subsection we will show how to obtain the results for the NMHV amplitudes using our 
generic formula for Tp^ given in eq. (|5.31| ). A coefficient of an Lj function will be of the 
form 

cT'^ = ^Pl,Pr ■ (6-6) 



"n,3 ~ 

Pl,Pr 



If the amplitude has three adjacent negative helicities they can only be positioned in two 
intrinsically different ways. In the first case, the massless leg is 2~ and the massive legs are 
r + 1+, ■ ■ ■ , 1~ and 3~, 4+, ■ ■ -r"^, i.e. the coefficients c^'g . In the second case, the massless 
leg has positive helicity and one massive leg contains two negative helicities and the other 
one, I. e. the coefficients c„ [ and c„ 2 • 

Consider the first case. In this case the paths are degenerate being from ai = 3 , pi = r 
to «! = 3, p2 = 4 for the right corner and from /3i = 1, cxi = r + 1 to /3i = 1 to (T2 = n on the 
left. 



'^0-2 = n 
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ai = 3, pi = r = 1, 0"! = r + 1 

right corner left corner 



The sum in eq. ( |6.6| ) then just becomes a single term, which can be computed along the lines 
of the previous section. The explicit expressions for c^'g are given above in eq. ( p.2| ) or in 
eq. ( |6.4[ ) for the two choices for P and P , depending on whether a scalar or A/" = 1 multiplet 
circulate in the loop. 

Let us consider now the second case. Both coefficients c„ 3 as well related to 



^m— 3,n— 2 



by parity and relabelling. We consider the coefficient d^'^ in more detail. The left 

. Am — 1)+) with m — 4 



massive cluster of legs of the corresponding triangle is (2 ,3 , 4+, . 
positive helicities. This coefficient is related by conjugation and relabelling to 



— m.l 



C, 



m— 3,n— 2 I 



n,n-3 l(123..n)-»{456...3),parity-flip " 



(6.7) 



The coefficient ^ can be computed along the lines of the previous section and gives two 



kinds of contributions. For both the right corner contributes with ao = 
and pi = n — 2 , which corresponds to a single path as displayed below. 



m — 2,ai = n — 3 




P2 = l 




^2 



n 




aQ = m — 2,pi=n — 2 
right corner 



l3o = Tn — 4,ai = n — 1 
left corner, k = 



0"2,3 



n 



/3i=m — 4,(Ti=n — 1 
left corner, k = 1 



For the left corner we have two cases which have to be distinguished. The first one comes 
from j3o = m — A,i3i=j,l32 = l,cri = n — 1 and a2 = n . One finds a term for the values of 
j = l,...,m — 5. The other contribution comes from the single argument Pi = m — 4 , j32 = 1 , 



ai = n — 1 and a2 = n . The sum over all paths that contribute to c^'^ is thus 



—m,l 
Cri,3 

-m,2 
^n,3 



^5'n,m— 1 ) 
~^^n,m— 1 ; 



(6.8) 



where gn,r and gn,r are given explicitly in eq. (|6.2|) or eq. ( |6.4|) for the two values of P and P 
corresponding to an A/" = 1 multiplet or scalar in the loop. Notice that the sum over j in the 
expression for Qn^r represents the contributions for the two cases k = and k = 1 together. 
It is striking that the computation of these coefficients is no more difficult than that of a 
tree-level calculation. 



7. Conclusion 

In the endeavour to compute the perturbative interactions within gauge theories, techniques 
which make direct use of previously computed amplitudes to generate new ones are usually 
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to be preferred. The loop-level unitarity method [38,39] or tree-level on-shell recursion rela- 
tions [11, 12] are examples of this. In this paper we developed an on-shell recursive method 
for obtaining coefficients of integral functions from previously calculated ones. We presented 
sufficiency criteria for avoiding spurious singularities and boundary terms in the recursion. 
These criteria apply to a large class of coefficients. 

We illustrated this with a non-trivial example, obtaining the complete n-gluon scattering 
with an A/" = 1 chiral multiplet circulating in the loop and where the external gluon helicities 
are "split", An{l~ , 2~, . . . {j + l)"*", . . . , n+) . For this same helicity configuration with 
a scalar circulating in the loop we also obtained all logarithmic terms in the amplitudes. 
Our starting point in the recursion were the cut containing parts of the five-point MHV 
amplitudes obtained in ref. [46]. With the supersymmetric decomposition, the contributions 
we have computed are pieces of QCD amplitudes. 

A particularly simple subset of the amplitudes obtained here is the NMHV case with three 
colour- adjacent negative hehcities and the rest positive. In this case, the M — 1 amplitudes 
agree with previous computations [21,29], providing a non-trivial check on our methods. The 
cut parts of the scalar loop contributions are new. Since the NMHV n-point amplitudes with 
an A/" = 4 multiplet circulating in the loop are also known [25], the only missing pieces for 
complete QCD amplitudes are the rational function terms. 

It would be very interesting to apply the ideas discussed in this paper to obtain other 
coefficients of integral functions appearing in one-loop QCD amplitudes. A new means of 
dealing with the rational function terms of loop amplitudes has also been given recently [2]. 
We may look forward to many new multi-parton one-loop amplitudes for use in collider 
physics. 
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